Propagation of surface plasmons through planar interface 
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We analyze the scattering of the surface plasmon incident at a planar interface between two 
dielectrics. By using the scattering matrix technique, developed by Oulton et al. [Phys. Rev. B 76, 
035408 (2007)], we calculate the transmission, reflection coefficients and radiative losses for oblique 
incident angles. We found that the transmission of a surface wave through a single interface between 
two dielectrics may be accompanied with radiation losses of 10-40 per cent of the plasmon energy. 



I. INTRODUCTION 



The excitation of a surface electromagnetic wave at the metallic interface - surface plasmonii^ - opens new ways in 
<^ i- nanophotonics^"^ and metamaterial physic^. One of the main constrain in using of surface plasmons is their short 
lifetime. It is well known^ii"— lii that a significant part of the plasmon energy is radiated when the plasmon is scattered 
at the surface impurity. A detailed quantitative analysis of the process of scattering and estimation of radiation losses 
is therefore important for understanding of propagation of surface plasmons. 
—I The most simple scattering problem is the transmission and reflection of a surface plasmon at the planar permittivity 

Y\ • step, created when the metallic surface is covered by two different dielectrics In the most simple scattering 
experiment, the metallic surface lies in the z = plane, and two dielectrics fill the z > half-space. The interface 
^ • between two dielectrics is given by a; = plane, so that the dielectric permittivity is 



g : '''^[el x>Q. 

^ . In this paper, we study the propagation of the surface plasmon through the planar interface between two dielectrics 
J>^' which cover the metallic surface. A modified method of Oulton et ai,^^ is used for the calculation of the transmission 
, and refiection coefficient and analysis of the radiation losses accompanying the plasmon scattering. Our data confirm 
that significant part of the surface plasmon energy is radiated in the process of single scattering, and energy losses 
increase when the angle of incidence increases. 
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\ II. SURFACE PLASMON AT THE METAL - DIELECTRIC INTERFACE 

o 

The surface plasmon propagates along the metal dielectric interface located in the z = plane. The intensity of 
the electric and magnetic field decays exponentially on both sides of the interface: h cx e"'^''^ for z > (dielectric) 
and h oc e"*'""*^ for z < (metal). The parameters Kd and k„i are given by the dispersion relations^^ 
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(sm is the metallic permittivity), and 

>< : ''I'^l ^ i^' ' ^ ° (3) 



Here, k\\ = y A:^ + ky is the projection of the wave vector into the xy plane, /cq — oj/c and c is the hght velocity. From 
Eq. (|2l3p we find explicit expressions for the components of the wave vector, 
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These equations, together with the Drude expression for the metallic permittivity, = 1 — w^/oj^, determines 
completely the frequency dependence of the wave vector of the surface plasmon. 

The surface plasmon is TM polarized, with magnetic field parallel to the metal-dielectric interface. The intensity 
of magnetic and electric field is of the form 

= AAo(-sin^,cos0,O) e*(fe^-+'=»?^-"*) x \ ^ J JJ^ (5) 
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FIG. 1: Relation between two angles 6t and 9a for various frequencies u of the surface plasmon and permittivity step et/ea ~ 5. 
In the case of the transmission a ^ b, the reflection angle 6b is always larger than that for the plane wave. Consequently, in 
the case of the transmission 6 — > a the critical angle 6cP is always smaller than the critical angle 9c for plane waves. 



and 

°ko \ (-iK„iCOs6',-m„sin0, -fc||)e+'''"^/£„i z<0. ^' 

Here, zq = yTio/eo is the vacuum impedance, and 9 determines the direction of propagation in the xy plane: 
cos 9 = kx/k\f^, sin0 = ky/k^^. Normalization constant Aq is specified in Appendix 1X1 



A. Snell's law for surface plasmon 



In the scattering experiment, we consider the metal covered by two different dielectrics, a and b, with permittivities 
£a and Sb- The interface between dielectrics lies in the yz plane x — 0. From the continuity of the y— component of 
the wave vector, kya = fcyj,, we find the relation between the incident and the refractive angle, 

sin 9b _ fe£ j Sb + £m 
sin 9a \l eb\l ea + 

Figure [1] shows 9b as a function of 9a for various values of the plasmon frequency w. The most important consequence 
of the relation ([7]) is the existence of the critical angle 9cp for the surface plasmon incident from the media with higher 
permittivity. No transmission of surface plasmon is possible when the incident angle 9b > 9cp. 



III. TRANSMISSION AND REFLECTION COEFFICIENTS FOR THE SURFACE PLASMON 



The transmission and reflection coefficients will be calculated from the requirement of the continuity of tangential 
components of both electric and magnetic fields at the interface x — Q between two dielectrics a and 6, 

E,a{x^{r) = E,b{x^Q+) and Hya{x ^ ir ) = Hyb{x (8) 

Since the intensity of the surface plasmon decreases exponentially in the z direction (Eqs. I5|6p . two plasmons on the 
opposite sides of the interface cannot satisfy the continuity relations ([8]) for all z. Therefore, we have to consider the 
full system of eigenwaves for the metal-dielectric interface. This system contains, besides the surface plasmon, an 
infinite number of plane waves. All plane wave have the same y— component of the wave vector, but differ in the z— 
component kz- In our numerical analysis, we use N plane waves with kza = (a/A^)fcmax and a = 1,2, ... ,7V. The 
upper cutoff fcmax is specified in Appendix |X| 

With the use of the plane waves, we express the continuity equations for the electric and magnetic field given by 
Eq. dH) in the for m 

JV JV 

^(^a — Aa)Ezaa = ^(^ct — Ba)Ezba (9) 
Q=0 Q=0 
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for the z-components of the electric field Ez, and 

N N 

+Af,)Hyap = Y,{Bp + Bi3)Hy,p (10) 
p=n 13=0 

for the y components of the magnetic field Hy. Fields Ez and Hy are given in the Appendix El Vectors Aa = A{kza) 
and Bp = B{kz/3), a, /3 = 0, 1, . . . , iV contain amplitudes of the surface plasmon (a, 13 = 0) and N plane waves in the 
left and right media, respectively. A, B {A, B) represent fields propagating to the right and to the left, respectively. 

All amplitudes can be calculated from the requirement that Eq. (|H]) must be fulfilled for any 2^'^. Another approach, 
which uses the coupling coefficients between electric and magnetic fields, suggested in Ref^ii, enables us to formulate 
the problem in terms of the 2(A^ + 1) x 2{N + 1) scattering matrix S 




(11) 



which relates the amplitudes of the incoming waves A and B with the outgoing waves A and B. Details of the 
calculation are given in Appendix [X] 

In the next Section, we analyze the case when the only incident wave is the surface plasmon propagating in media 
a. Then A^ = 5^^ and 5 = 0. From Eq. ([Tl]) we obtain the transmission and refiection coefficients for the surface 
plasmon, 

ra^b = |SgS|2 and Ra^a = \SZ?- (12) 

The components SJ^g and SJ^q determine radiation losses due to the scattering of the surface plasmon. Among all 
plane waves, only those with real kx radiate the energy in the x direction. Since fc^^.^ = k^Sa — ky — k'^zcn have 
that kaxa is real only for a smaller than certain integer Ua- Similarly, khxp is real only when /3 < n^. Total radiation 
losses are therefore obtained as Sa = Saa + Sba , where 

ria rih 

Saa = ^\SZ\^ and Ska^J2\S%\'. (13) 

a 

The conservation of the energy requires 

Ta^b + Ra^a+Sa^l. (14) 

Physical meaning of other components of the scattering matrix is obvious. For instance, the element Sq^ gives the 
amplitude of a surface plasmon, excited in the media a by a plane wave /3 incident to the interface from the media b. 



A. Normal incidence 



Figures [2] - [4] show scattering parameters of the surface plasmon for the case of normal incidence. In numerical 
calculations, we use fcmax given by Eq. (jA3p and number of plane waves varies between N = 100 and TV = 1577. 

The transmission and reflection coefficients as well as radiative losses are given in Fig. [5] for various values of the 
permittivity steps Sb/Sa- Our data agree with results of RefJ^ and confirm that a significant part of the plasmon 
energy is radiated by plane waves. Radiation losses increase when the permittivity step increases. On the other hand, 
scattering coefficients depend only weakly on the plasmon frequency (data not shown). 

As the test of numerical accuracy of the method, we used amplitudes A and B, obtained from the scattering matrix, 
and reconstruct the electric and magnetic fields on both sides of the x = plane for the permittivity step £b/£a = 5. 
Figure [3] confirm that the tangential components of Ez and Hy are indeed continuous at the interface. 

In Fig. |3]we present the amplitudes of radiated plane waves, |SJ^qP for the plasmon incident from the media a and 
h. The data confirm that the energy is mostly radiated in the direction of incoming plasmon. As shown in the right 
figure, the radiation possesses the sharp maximum in the direction of the critical angle for planar waves Oc- 
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FIG. 2: Transmission T , reflection R and radiative losses S for the normal incidence of the surface plasmon at the interface 
between two dielectrics. The surface plasmon is coming from the dielectric a (left) and from the dielectric 6 (right). 



B. Oblique incident angle 

Figures [5] and |6] show how the scattering coefficients depend on the incident angle. The plasmon is approaching the 
permittivity step either from the left or from the right side of the interface. 

For the scattering from the side with lower permittivity, ea < £b, all coefficients depend monotonously on the 
incident angle. More interesting is the case when the plasmon approaches the interface from the side with higher 
permittivity Sb- The transmission, T^^a decreases to zero when 9 — >■ 6cp, but the reflection R does not increase to the 
unity. We explain this behavior by the presence of "evanescent plasmon" in the media a. Although the x— component 
of the plasmon wave vector kax is imaginary, the intensity of the field on the left side of the interface is non-zero (even 
larger than for smaller incident angles). This field must be compensated by plane waves which radiate energy. 

For higher permittivity contrast Sb/Sa, we found that the reflection even decreases when the incident angle increases 
above the critical angle 9cP- This decrease is accompanied by higher radiation losses. As shown in Fig. [6] entire 
plasmon energy can be radiated when 9 > 9cp. Radiation losses have a maximum for the incident angle larger than 
the critical angle. The "total reflection" (i? = 1) takes place only for angles much larger than the critical angle for 
the surface plasmon. 




FIG. 3: Test of the continuity of electric and magnetic field along the interface a; = 0. Solid lines and symbols represent fields 
for X ^ 0~ and a; — >■ 0'*', respectively. Left figure show real part of both E and H, and two other figures show imaginary parts 
of fields. Dielectric permittivities are Ea ~ £b ~ 5. N — 1577 plane waves were used with with maximal z— component of 
wave vector given by Eq. (|A3|) . The maximal value of z is given as Zn-m^ = JtljKda- 
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FIG. 4: Angle distribution of radiative losses E{kz) (in arbitrary units) given by scattering matrix elements S^q for the 
scattering of the surface plasmon with the frequency uj = 0.23a;p scattered at the interface between two dielectrics with 
permittivities Sa = 1, £5 = 10. fcoa = koy/e^ and fcot = koy/ib- Left (right) figure shows radiation in the medium a (b), 
respectively. A and B determines the medium from which plasmon is coming. Sharp maximum in the right figure shows the 
scattering in the direction of the critical angle 9c for plane waves. 



IV. CONCLUSION 



We analyzed quantitatively the scattering of the surface plasmon at the planar interface between two dielectrics 
which cover the metallic surface. The transmission, reflection coefficients and radiative losses were calculated for the 
normal and oblique incident angle. We confirm that the radiation of plane waves causes significant scattering losses: 
for normal incidence, the transmission through the single interface might cost 20-40 % of the plasmon energy. The 
reduction of these losses represent the challenging problem for the theoretical research. One possible way how to avoid 
this problem is to cover the metallic surface by anisotropic metamaterial instead of a dielectric"'^'^. 

We analyzed how the transmission and reflection coefficients depend on the incident angle. While the transmission 
to the dielectrics with Eb > £a brings no surprising result, the transmission in the opposite direction exhibits non- 
monotonous dependence on the incident angle. The transmission coefficient decreases to zero when the incident angle 
increases to the critical angle OcP for the surface plasmon. However, the reflection does not reach unity for the critical 
angle, because the significant part of the energy is radiated. 

In the present analysis, we used real (losless) metallic permittivity. This is consistent with the formulation of the 
scattering experiment, in which the incident wave is coming from the infinity. Nevertheless, we verified that realistic 
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FIG. 5: Transmission T, reflection R and radiative losses S as a function of incidence angle for the frequency uj = 0.23ajp and 
dielectric interface Sa ~ 1 and et ~ 2. Left: scattering a ^ b. Right: scattering b ^ a. Vertical dashed lines show critical 
angles 9cP for the surface plasmon and 9c for plane waves. 
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FIG. 6: Transmission Tt^a, reflection Rb^b and radiative losses Sb as a function of incidence angle for the frequency uj = 0.23ajp 
and interface Sa ~ i and £5 = 6 (left) and Sb = 10 (right). Vertical dashed lines show critical angles 6cP for the surface plasmon 
and 6c for plane waves. The transmission is zero when > Ocp. The reflection decreases for 9 > 6cP in favor of the radiation 
losses Sb- Total reflection is observable only for angles 9 ^ 6c- 



losses, given by small imaginary part of the metallic permittivity in the Drude formula, do not influences the scattering 
coefficients. 



Appendix A: The method 

Since the single surface plasmon cannot satisfy the continuity relations ([8]), along the dielectric interface, a complete 
set of plane waves must be included into the scattering procedure. Each plane wave is given by a superposition of 
the wave incident to (cx e"'*^"^^^), reflected from (cx re"'"*'^''^^) and transmitted through the metal-dielectric interface 
(cx ie"*'^™^^). Here kdz and kmz are the z— components of the wave vector in dielectric and metal, respectively. The 
reflection amplitude r for the metal-dielectric interface is given by 

Sdkmz ^mkdz (Al) 



dz 

and t = \ — r- In numerical calculation, we consider N plane waves with different values of the z-component of the 
wave vector 

fcd.a = -^xa, a = 1,2,..., TV. (A2) 



Here, fcmax is the largest allowed value of kz- We choose 



where Sd = min(£a,£6)- This choice guarantees that all waves transmitted from the dielectric to the metal decrease 
exponentially, so that no plane wave propagates inside the metal and the z— component of the wave vector in the 
metal. 



max? 

(A4) 
is imaginary. 

In what follows we need explicit form of the the plane wave for the interface metal-dielectric a the fields Ez and 
Hy. Neglecting the phase factor exp[i(fcj.a; -\- kyU — ujt)], we have 



and 
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with kmzaa given by Eq. (jA4p . Similar expression can be written for the interface metal-dielectric b. 

To solve the system (|9ll0p . we multiply both sides of Eq. ^ by Hyaa and Eq. ([TU]) by E^bp and integrate over z. 
Since all waves possess the same the y— component of the wave vector, the coupling coefficient between two plane 
waves is given by the integraliiii^ 



Cab 
a/3 



dzEznnH, 



zaa^yb(5 : 



(A7) 



where ^ is infinity when a ^ b and ^ = Zmax = {2tt /ki^a.^)N for a — b. The choice of the spatial cutoff Zmax enables 
us to express all integrals (|A7|) in terms of dimensionless parameters. We find that diagonal elements of the matrix 
C reads 



Cab 
aa 



-MaaMbaTTT. Z i^ic + na) ■ 



The requirement C^^ = Sa/s determines the norm J\fa. 



J^a, 



zq 2ttN 

^aa '^a 



The off-diagonal elements read 



— —iNaaJ^bjskbxp- 



ko kb\\ £m(fcm2f,^ ^mzaa) 

Similarly, diagonal elements for two plasmons. 



kmzbl3 k^^da (^mzfc/3 kmzaoi) 

^b 



Cab 
nn 



dzBzahyb = A/'aoA/'bOT^T^/Cha; 

ko kb\\ 



1 1 



1 



£a l^ad + l^bd £m 1^, 



am f^brn 



determines the normalization constant for the surface plasmon, 

A/'aO = 



I kg K da "^Sa^r, 



Zokax £rn ^a 



(A8) 



(A9) 



(a,^>0). (AlO) 



(All) 



(A12) 



Using the form of the electric and magnetic field of the surface plasmon, we obtain the coupling coefficients between 
the surface plasmon and the plane wave in the form 



ah 



+ 00 



dzCzaliybP = ■N'aO-N'bp-r- T^^bxli 
ko Kb\ 



1 



1 



Had - ikbzfS Had + ikbzf) 



1 np - I 



bp 



(A13) 



Finally, we obtain two sets of -I- 1 linear equations for unknown amplitudes A, B, A and B: 

A-A_= C^{B-B) 
B + B = C (A + A) 

which can be rearranged into the form 



with 



S''^ = -(l + CC^)-i(l-C 



/gbb 


o6a \ 
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gba 


= (1- 




gaa 


= (1- 
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(A15) 



(A16) 
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